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Abstract: We consider a collapsing relativistic spherical shell for a free quantum 
field. Once the center of the wavefunction of the shell passes a certain radius r^, the 
degrees of freedom inside are traced over. We show that an observer outside this 
region will determine that the evolution of the system is nonunitary. We argue that 
this phenomenon is generic to entangled systems, and discuss a possible relation to 
black hole physics. 
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1. Introduction 

Consider a free massless scalar field in a pure state describing a spherically 
symmetric inward collapsing shell. Imagine that once the shell reaches a certain 
radius r^, an observer located at r > has no access to the region r < r^. This 
implies that the state seen by this outside observer is a mixed state described by the 
density matrix pout = Trr<rjV') (V'l- We shall investigate the time evolution of the 
state Pout- 

There are general arguments, which we discuss below, that show that a generic 
entangled state will seem to evolve in a nonunitary way if part of the Hilbert space 
is inaccessible — a property that leads to dissipation effects in many body systems. 
In this paper, we shall show that the eigenvalues of pout are time dependent. This 
implies that pout does not evolve in a unitary manner (see figure |l|). To see this, 
consider a quantum system prepared in an initial state defined by the density matrix 
p(0). If the evolution is unitary, then the density matrix p{t) describing the system 
at a later time is related to p(0) by pit) = W (t) p{0)U (t) . Since this is a similarity 
transformation, the eigenvalues of p{t) must be the same as those of p(0), and are 
consequently time- independent. Conversely, consider a density matrix p{t) whose 
eigenvalues are time dependent. In this case, for the same reason, a unitary operator 
U(t) such that p(t) = W(t)p{0)U(t) does not exist. Therefore, the time evolution of 
a system will be unitary if and only if the eigenvalues of p are time independent. 

Previous investigations of the relationship between entanglement, entropy and 
area |jl| have shown that the von Neumann entropy S = Tir^rsiPout^^ Pout) of an 
initial vacuum state |0) is proportional to the surface area of the region r < r^, 
suggesting a connection with black hole entropy. This observation (see also 0) has 
been studied in the literature from various perspectives ^, ^, |^, |^, ^ |TU|, |TT|, 
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131 , [T^, |T5|, |T6|, ^7} |T8|, |T9[ (see ^] for a review). In flat space, it can be shown that 
the entropy (and other thermodynamic quantities) will scale as the surface area of the 
inaccessible volume even when it is nonspherical ||I2|, 115], 1131. This area dependence 



may be related to thermodynamic properties of Unruh radiation |22, 23 



The flat space arguments relating entanglement entropy and area may be ex- 
tended to a black hole background [El, |2S|, BH] . Since this work is key in understand- 



ing the relation between our results to black hole evolution, we shall briefly review 
it. Consider, for simplicity, an eternal uncharged and non rotating black hole. This 
black hole has two asymptotic regions separated by a horizon. As a result of this 
causal structure, an observer in one asymptotic region has no access to the degrees 
of freedom in the other asymptotic region. In contrast, the global vacuum of this 
system (the Hartle-Hawking vacuum) extends to both asymptotic regions. An ob- 
server in one asymptotic region will observe the state obtained from the vacuum after 
degrees of freedom in the other asymptotic region are traced over — the degrees of 
freedom "behind" the horizon. The resultant state is not a pure state but rather a 
mixed state which turns out to be precisely the density matrix seen by an observer 
restricted to one asymptotic region at the Hawking temperature. As in the flat space 
case, the entropy associated with this state is proportional to the surface area of the 
horizon. This mechanism, which is a result of the entangled nature of the vacuum 
state, allows one to interpret the black hole entropy as coming from entanglement, 
an interpretation which is consistent with the string theory evaluation of black hole 

for a recent discussion 



,29 



see also |31, 32, 33, 34 



entropy |§ 

of black hole entropy and entanglement in the context of AdS/CFT.) 

Thus, following the current work is a natural next step in relating features 
of entanglement to black hole physics. While it is generally believed that black holes 
have entropy proportional to their surface area ||38| there is, at the moment, some 
debate regarding their time evolution. 



In it was argued that the evolution in time of a black hole is nonunitary (see 



41| , ^ |45| , |49| , |50| for recent discussions and implications). 



This is a troubling feature of black holes since it clashes with our understanding of 
quantum mechanics. For instance, if the evolution is nonunitary a pure state (for 
which the density matrix has a single nontrivial eigenvalue) may evolve into a mixed 
state. From the entanglement point of view this behavior is not in contradiction 
with the principles of quantum mechanics or any laws of physics, rather it is a simple 
consequence of the entangled nature of the system. If one does not have access to 
part of the Hilbert space then probability flowing into the inaccessible region will 
seem to disappear leading to nonunitary evolution. 



2. General arguments 

There are several approaches one may take when trying to exhibit the nonunitary 



-2- 



evolution of an entangled system. Many of these have been studied in the context 



of dissipation of energy from a system to a heat bath (see for example and 
references therein) and may be directly applied to our discussion. We review them 
here for completeness. 

We would like to exhibit the nonunitary evolution of an initial pure state p(0) = 
!'?/')('?/' I seen by an observer who has no access to a certain region of space. We 
evolve our initial state in time using the Liouvillian operator L = —i[H, ], so that 
p{t)=Cp{t). 

We shall split our Hilbert space of states into two parts, which we will refer to as 
the 'in' region and the 'out' region, such that the Hamiltonian can be written as H = 
Hout + Hi + Hin with Hi an interaction term. Since C is linear, this correspondingly 
results in £ = Cout + Ci + Cin- 

Consider an observer restricted to, say, the 'out' region of space, and let R 
be an operator that restricts the density matrix to this region. This corresponds 
to "integrating out" the unobservable degrees of freedom. (In the path integral 
formalism, R is given by eq. ( |3.2| ) below.) We assume that the 'in' and 'out' regions 
are fixed, so that R is time-independent (though a generalization to a time dependent 
R is straightforward). 

Suppose that at some time t = the state of the observer is given by Pouti,^) = 
Rp{0). The state pout will then evolve in time according to Poutif) = Rp{t). Note 
that, once we restrict observables to the 'out' region, a further restriction will not 
change the state of the system, therefore R^ = R. 

To find Pout, we act, in turn, with (/ — R) and R on p = Cpout + '^{I — R)p, 
obtaining equations for pout and {I — R)p. Plugging the (formal) solution to {I — R)p 
into the equation for pout "we obtain the Nakajima-Zwanzig [^, [ 
master equation 



54 form of the 



Pout{t) = {RC)pout{t) + / {RC) Lexp 



{{I-R)C) dt' 



t-t 



+ (RC) lexp (^j^ (/ - R)Cdt"^ (/ - R)p{0). 



{{I - R)C) pout\,_fdi+ 



[2.1) 



To simplify (|2.1|) we note that the restriction R commutes with time evolution in 
the 'in' and 'out' regions [R, Cout] = [R, Cin] = 0. Also, we assume that, CinR = 0. 



Finally, since i?exp 



((/ - R)C) dt' ((/ - R)Ci) (t - f) = 0, we find 



Pout{t) =R{Cout + Cl)pout(t) + 

+ [ {RCi)lexp [ {{I-R)C)dt' 
Jo \_Jt-t 



{{I - R)Ci) poutl_^dt+ 



+ RCiit) exp (^j^ (/ - R)Cdt"^ (/ - R)p{0). 



[2.2) 
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The first term in eq. ( p.2|) is the standard term that one expects for unitary time 
evolution of the state seen by the observer in the 'out' region. The second term 
represents the (time retarded) oscillations flowing into the 'out' region from the 'in' 
region that is being traced over ("beyond the horizon"). The last term represents 
contributions leaking out of the inaccessible region from those initially present. In 
most many-body applications, it is arranged or assumed that (/ — -R)p(O) ~ 0, so 
that this term can be neglected. 

The interaction term in the Hamiltonian, Hj, is essential for generating the 
nonunitary evolution. In a field theory setting, this interaction term includes the 
spatial derivative coupling fields across the boundary. Thus, it exists even for a free 
field theory. 



3. Explicit construction of nonunitary evolution. 

We wish to explicitly exhibit the nonunitary evolution measured by an observer who 
does not have access to part of a system. Consider the time evolution of an inward 
collapsing spherical shell of a free massless field in 3 + 1 dimensions from the point of 
view of an observer who does not have access to the region r < r^. We choose such 
a configuration for two reasons: (a) to make suggestive contact with gravitational 
collapse, and (b) since any choice of an initial state that is an eigenstate of the 
Hamiltonian (such as the vacuum state) will not exhibit any time evolution due to 
its stationary nature. 

Let \E'fc[0(a;)] be the wavefunctional of a single particle with energy i^^ = of a 
free massless scalar field (t){x) in three space dimensions. An arbitrary configuration 
^ for such a particle is specified by = Ylik fk^k- We construct a collapsing shell 
by choosing {fk] such that only s-wave modes are excited and such that the initial 
position space wavefunction is centered at some radius rj and is collapsing toward 
the origin. To obtain explicitly, consider the wavefunction 

\^) = / /(r)0(x)|O)rf3x (3.1) 



where A is a normalization constant determined by the condition = 1 and /(r) 

is real. The average momentum of the shell, p, is set by substituting /(r) e*^^/(r). 
The set {fk} may be obtained from the expansion of the field in terms of ladder 
operators. 

The density matrix of the region seen by an observer in the exterior is given by 



Pout{4>out,(j)'out) = j D(j)in'^{(j)in,4>out)'^{(j)in,(j)'outy- (3.2) 

A similar density matrix was constructed in a somewhat different context in EH . 
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As in [[T], 0, one can calculate the functional integral ( |3.2|) by discretizing space 
on a radial lattice. We expand the scalar field in partial wave components 



4>i,m{r)=r / Zim{0,^)(f){x)dn, T^i,m{r) = r / Zim{9,v)i^{x)dVL 



(3.3) 



where r = |x| and Zq^ = Fq.o, Zi^m = \/2Rey/,m for m > 0, and Zi ^a = v^ImF/^m for 
m < 0. Zi m are orthonormal and complete. This allows us to write the Hamiltonian 
as = ^ Hi m with 



l,m 



2 Jo 



dr 



(3.4) 



We discretize the radial coordinate on a lattice of spacing a and of size L = 
{N + l)a with boundary conditions such that the fields vanish at r = L. The 
expressions for the Hi^^s ( |3.4| ) become those of Hamiltonians of coupled harmonic 
oscillators. 



Hl 



1 ^ 

-y 



T^lm,i + (i + nf 



Pl,m,j+1 
J + 1 



+ 



la 



(3.5) 



Note that all the oscillators are coupled due to the spatial derivatives of the field. 
This coupling generates the interaction Hamiltonian described earlier. 

The wavefunction for a single particle shell is constructed as in the continuum 

m 

2 ^ 

l^) = :^E^/^'^OA.|0). (3.6) 

We wish to write this in a single particle eigenbasis of the Hamiltonian {aj|0)}^]^. 
Let Hi = 5'°-0o,oj) where 5"° diagonalizes Kq q defined in (|3.5| ). We find that 



A- 



1 ^ 



|0) 



(3.7) 



with 



3,f 



J J 



In ( |3.71 ) we have defined {VlPY = S^Kq^qS'^^ . As in the continuum, one can generate 
a collapsing shell by substituting fj — > e^^^ fj. The time dependence of the state 
can be read off from \il),t) = e~'^^^\il)) 
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The wavefunctional on this discretized space is given by 

with No = y = S(f) where S diagonahzes K, (K)" = and Hi{x) = 2x 

is a Hermite Polynomial. In what follows it will be useful to compare our results to 
those obtained by studying the ground state wavefunctional 

^o({0^,n^J}) = iVoe-^^^^. (3.8) 

We shall consider a general linear superposition of single excited oscillators 

^ = $^/fe^fc (3.9) 

k 

= ^(^205^fi^/)e-^'^^^^. (3.10) 

For convenience, we define 

vk = (^fn-^s^^. (3.11) 

Note that f„ = ^j3fj{^~^^^)j,n- This allows us to write pout = / '4''4'*d(j)in as 

Pout{w,w') = 2 !{v (t>W ■ c/'OliVolVi'^^^ V^^'^^^'d";^, (3.12) 



where we have used = (z, vo) , and </>' = (2;, w') . 

The right-hand-side of eq. (|3.12 ) is a Gaussian integral over the first n coordi- 



nates. Using the notation 

5^ C 



^^=fdf), (3.13) 



V2 



we find 

Pout{w,w') =2po^out{w,w')x 



- {v2 ■ w)-vlA~^B{w + w') - ■ w')-viA'^B{w + w') + 
2 2 



+ iv2 ■ w){v* ■ w') , (3.14) 
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where 



P0,out[W,W 



1T2 
A 2 



^{w + w')B^A~^B{w + w') - ]^{wCw + w'Cw')^ 

(3.15) 

One can check that 



see also 111 for an 



is the density matrix corresponding to the vacuum state 

Trpo.out = T^^pout = 1- 

The eigenvalues of po^out may be obtained numerically [|l| 
analytic approach). An approximation for the eigenvalues of the first excited state 
has been discussed in 0. Since we are interested in exhibiting the nonunitary time 
evolution of the state measured by an external observer at r > r^, we shall evaluate 
Tr(p^„i); if Tr(p^„J is time dependent then, since Ti^pout) = 1, at least two of its 
eigenvalues are time dependent. 

To calculate Tr(p^„J we carry out another Gaussian integral. For the ground 
state we find 

\K\ 

(3.16) 



Tr(pL«J 



where we have defined 



T 



C 



IB'^A-'B 

1 UT A~lj^ 



-{B^A 



\A\\T\ 



-\B^A-^B 
C-lB^A-'B 



(3.17) 



If K were block diagonal such that B = 0, then we would have had |T| = \C\ 



so that Tr(pQQ^J = 1. We expect this since B 
oscillators are not entangled. 
For our collapsing shell (|3 



implies that the 'in' and 'out' 
we find, after some algebra, that 



Tr(p 



where 



2 1 

out) 



A-^ 




u 



B 



f + \vUBvf 



Tr(p 



2 1 

O.out) 



Ua 



'\A-^B{oA + ob)B^A~^ -].a-^b{oa + Ob] 
\ -\{oa + ob)B^A-^ oa,b 

and OA and ob are defined through 



(3.18) 



(3.19) 



OA Ob 

Ob oa 



(3.20) 



(where T is given in (p.l7|) , and the block matrices A, B and C are given in ( p.l3|) ). 
As in the ground state case, one may check that Tr(p^^J = 1 if i? = 0. 

In principle, the time dependence of pout can be extracted from ( |3.18|) . It is clear 
that if / (recall v = fQ^/^S) has a single entry (implying that the state of the system 
is an energy eigenstate) then Tt{pI^^) will be time independent. The same applies 
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to the ground state Tr(pQQ^J (|3.16|) . In fact, as discussed earlier, any stationary 
state will generate a time independent pout-, meaning that the time evolution will be 
unitary. (Curiously, the same applies to coherent states (see 0).) To obtain the 
time dependence of Tr(p^^J, we shall resort to numerical methods. 

In figure 1^, we have plotted the radial wavefunction J2j J fji^~^^^)j,n as a function 
of the lattice location, n, and time, t, for the time interval [4a, 15a]. We used a radial 
lattice of size + 1 = 50. The initial wavefunction was chosen such that at t = 
it is a Gaussian shell localized at = 33a with width a = -y/Sa, and momentum 
—la. The thick (red) line specifies the location of the maximum of the wavefunction. 
Once the center of the wavefunction passes = 25a (at t = 8.84a) we trace over the 
region r < r^. In the figure, this region is bounded by a semi-transparent membrane. 
The value of Ty{pI^^) at t > 8.84a is also plotted. 

As the wavefunction enters the inaccessible region r < 25a, Tr(p^^J starts in- 
creasing. This increase is to be expected: once the wavefunction is 'almost' com- 
pletely inside the region r < r^, then the one particle excitation of the density matrix 
Pin would describe an 'almost' pure state implying Tr(pf„)/Tr(pQ j„) = 1. Since pi„ 
and Pout have equal eigenvalues (up to zeros), we get the behavior depicted in the 
plot. 

If the tracing procedure is carried out before t = 8.84a, then one finds that 
Ti'(pLJ/Tr(Poo«t) drops smoothly from ~ 1 (when the wavefunction is localized 
at r > Ts) to ~ 0.4 at t ~ 8.8a. Similarly, one can extend the analysis to t > 15a, 
whence the wavefunction "bounces back" from the origin: as the wavefunction leaves 
the r > Tg region, another sharp change in Tr(pQ„J/Tr(pQ is noted. At very long 
time scales, the wavefunction will bounce back and forth from the origin and IR 
boundary of space. In this case one will find a semi-periodic behavior for Tr(p^^J. 
In general, we observe that Tt{pI^^) = a;(t)Tr(pQ (see eq. ( p.l^ )), where the 
coefficient a approaches unity when the wavefunction is localized either in the r < 
region or when it is localized in the r > region. 

We point out that since the wavefunctions corresponding to 0o,o,i|O) are orthog- 
onal under the Klein Gordon norm then they are not orthogonal under the Dirac 
norm. Therefore, the plot only gives a suggestive description of the region where 
the wavefunction is localized. Also, since the scalar field is massless, a zero momen- 
tum s-wave will split into an outgoing shell and an ingoing shell. The amplitude of 
the outgoing shell is small initially because the wavefunction has been given a large 
initial ingoing momentum. 

4. Discussion 

The main point of this exercise is to emphasize that nonunitary time evolution of a 
density matrix is an ubiquitous phenomenon and does not signify a breakdown either 
of unitarity or of quantum mechanics. Whenever an observer has a limited domain 
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lattice location [a] 



time [a] 



Figure 1: A plot of the radial component of a collapsing spherical wavefunction as a 
function of the lattice location, n, and time, t. The radial lattice has 50 points of spacing 
a. At t = the wavefunction is a Gaussian shell localized at 33a with width a = ^/ba, 
and momentum — 1/a. The thick (red) line specifies the location of the maximum of the 
wavefunction. Once the center of the wavefunction passes = 25a (at t = 8.84a) we trace 
over the region r < (bounded by a semi-transparent membrane) . rJ^^^°"*\ is also plotted 
for t > 8.84a. 



of observation, she will describe observables in terms of a density matrix resulting 
from tracing over states outside that domain [56]. If her states are entangled with 
states beyond her domain, then she will experience a density matrix corresponding 
to a mixed state. It may be static, or it may be time-dependent. If time-dependent, 
it may or may not be describable in terms of the dynamical variables associated 
with her limited domain of observation, depending on the situation. Although the 
global state may be simple and involve unitary evolution, the time dependence of 
the density matrix of a subdomain can be quite complicated and appear nonunitary. 

What implications might this have for gravitational collapse of matter to a black 
hole? Once a horizon forms, all observables can be described in terms of a density 
matrix associated with the states within a causal volume since, obviously, states 
beyond the horizon are unobservable. This situation is inherently time-dependent 
due to black hole radiation. By energy conservation, the horizon should shrink as 
this radiation is emitted. In our view, the time dependence of the density matrix 



describing the region outside the horizon is very hkely to appear to be nonunitary. 
Globally, we expect that the evolution is unitary and, if the matter started in a pure 
state, the total entropy will remain zero. In this case, the apparent entropy of the 
black hole would be entirely the result of entanglement. 

One major difference between our setup and a collapsing black hole is that our 
example is not singular and hence does not contain an analog of what happens to 
matter as it falls into a region of high curvature, where a classical singularity appears. 
Yet , if matter does not disappear from the universe that remains when the black hole 
has completely evaporated, one may speculate about the time evolution of the black 
hole based on this example. On one hand, one will find that from the point of 
view of a causal observer the evolution of the black hole will seem nonunitary. On 
the other hand, from a global viewpoint the wavefunction will evolve in a perfectly 
unitary manner. The dichotomy between the two observers is only a result of one 
of them restricting her observations to part of the Hilbert space, much like the case 
of a system coupled to a heat bath. At the end of the day, when the black hole 
has evaporated, there should be a unitary transformation from the state before the 
horizon formed to the state after the horizon has disappeared. 
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